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Abstract 

It is shown by a straightforward argument that the Hamiltonian generating 
the time evolution of the Dirac wave function in relativistic quantum mechanics is 
not hermitian with respect to the covariantly defined inner product whenever the 
background metric is time dependent. An alternative, hermitian, Hamiltonian is 
found and is shown to be directly related to the canonical field Hamiltonian used in 
quantum field theory. 



1 Hermiticity of the Dirac Hamiltonian 

We start with the Dirac Lagrangian in a gravitational background (see P) 

L = ^(V^f D.V' - d.V'tV - ^V'V'), (1) 

with 'Diij) = diip — IF'^^'jCTabV' cinci 7* = 6^7", where 7'* are the constant Dirac matrices and 
aab the Lorentz generators. Our notations are identical to those used in j.2 . Especially, 
we use latin indices from the begin of the alphabet (a, 6, c . . .) for tangent space quantities 
and from the middle k . . .) for curved spacetime quantities. The Lorentz connection 
-pafe^ (Yab^ _ — r^"J and the tetrad transform under a local Lorentz transformation e"'' 
as Sr^-^i = —die"''' — T"■^^e'''' — ^''^^e"■'', Se^ = e^e^. The Lagrangian (1) is invariant under 
such transformations if ijj undergoes the gauge transformation ip — > exp [{—i/4)e"'''aab]4'- 
The Dirac equation is derived in the form 

i'^^'Vi'ip = mip, (2) 
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where Vj'?/' = diifj — ^r°-'^^aabip- The connection F^^'j is the connection that effectively 
couples to the spinor field. If we work in the framework of general relativity, then r^''^ = 
r"^j(e), i.e., the connection is not an independent field and can be expressed in terms of 
the tetrad and its derivatives. In that case, F"''- = T"-'^- (and Dj = Vj). If the connection 
is considered to be an independent field (Poincare gauge theory), then we have T"-^- = 
r"'^^{e)+K"-\, where K"'\ is the contortion tensor. However, only the totally antisymmetric 
part of K"'''^ remains in the field equations. Indeed, it is not hard to see that the remaining 
parts cancel already in the Lagrangian (1). Thus, the effective connection is given by 
r^^j = T°'''-{e) + K"''^i, where K"'\ = eciK^""^^^- Those differences between general relativity 
and Poincare gauge theory are completely unrelated to the discussion of this paper, and 
the use of the symbols Vj and F^^j is a convenient way of treating both cases at the same 
time. 

For the study of the time evolution of ip^ the Schroedinger form of the Dirac equation, 

= idt^, (3) 

is especially adapted. We use greek indices a, /3 . . . to denote the spatial part of a,b . . . 
(i.e., a = (0, a)), and indices /x, z/ ... for the spatial part of the spacetime indices i,j,k... 
and t for the time component, (e.g., m = (t, /i)). With the metric gik = e°;elriab and using 
7*7* = (yf**, we find explicitely 

H = 4(7*^ - - Ir^Wb. (4) 

In the fiat limit (F"*^ = 0, = 5"), this reduces to the well known expression H = 
7°m — i'~i^'~i^d^ = [3m + a ■ p. 

The Dirac Hamiltonian in curved spacetime has been used in the past, see, e.g., OS 
inilSlIZ]; with the focus mainly on the coupling to gravity in the weak field limit. Although 
it has been shown that H is not generally hermitian j^l E] , this is not widely known, since 
in most relevant cases, the gravitational fields are static, and H is hermitian. In this 
article, we will analyze in detail the hermiticity properties of H. 

Before anything can be said about the hermiticity, we have to define a scalar product. 
One could try to construct this product such that certain conditions, hke covariance, 
or hermiticity of H are a result of the definition. However, there exists an even more 
physical requirement that leads straightforwardly to a consistent scalar product, namely 
the possibility of the normalization of the wave function. 

Indeed, using the Dirac equation (2) and its hermitian conjugate [hermitian conjugate 
refers here only to Dirac space, i.e., the transposed complex conjugate equation), one 
derives the following conservation law 

9,(e^7V) = 0, (5) 
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where e = det = \J—9- Integrating and omitting boundary terms, we get 

= ^ y d^x^-Jy/^ d^x. (6) 

Therefore, obviously, the only way we can obtain a constant in time normalization is by 
using the inner product 

{^1,^2) = / 4(lW)^2 d'x, (7) 
or, in an explicitely covariant form 

{il^i,ilj2)= I ^ dSi. (8) 



This product is covariant and positive definite and the usual properties (linearity, (■^i, ■02) = 
(V'2!V'i)*) ^re easily verified. In the flat limit, where 7* = 7°, it reduces to the special 
relativistic expression (V'i,V'2) — / V'l d^x. There is nothing else we could possibly 
require from an inner product. 

With (7) or (8), it can be directly checked whether H from (6) is hermitian or not, i.e., 
whether we have {ip,Hip) — {Hip,ip) or not. This involves, however, quite a few partial 
integrations and a lot of 7— matrix shuffling. Fortunately, there is a much easier way to 
get to the answer. 

We start from equation (6) and write 

o^y" [dtilj^^W^ + ^p^^WdM + [{dt^/^)i^hVi^ + V^i^h''{da')i^]]d^x, 

then use 7*^ = 6*^7", as well as the Dirac equation Hip — idti) and its conjugate ip'^H^ ~ 
—idtip^. (We remind once again that W is the hermitian conjugate with respect to Dirac 
space only, but we take the additional convention that its differentiation operators act to 
the left.) After some simple manipulations, we finally find 

= iP) - {iP, Hil;)] + (V', [dt In v/=^ ]^), (9) 

showing that H is, in general, not hermitian. In view of g^* det Qik — detg^^jy, this will 
be the case whenever the determinant of the three dimensional space metric g^j^i, is time 
dependent. Moreover, we can define a hermitian operator 



H = H + -dt\n^/^\ (10) 



satisfying 



(11) 



Thus, it turns out that the Hamiltonian in a gravitational background is not generally 
hermitian. We should note, however, that in most cases, det (7^,^ will be time independent. 
This holds notably for static and stationary spacetimes, but also for gravitational waves 
when the usual gauge conditions are used in the linear approximation (see [8j). Coun- 
terexamples can be found, e.g., in exact wave solutions of general relativity, as described 
in [Hj, §109, and, most importantly, in cosmological solutions jS]. 

Let us also remark that, as we have claimed at the begin of the section, the final 
conclusion is completely independent of the eventual presence of torsion fields, since the 
non-hermitian term in (9) depends only on the spacetime metric. 



2 Field Hamiltonian 

The main problem we face now is the interpretation of both H and H. On one hand, 
we have a non-hermitian operator, with, in general, complex eigenvalues, but with a 
straightforward interpretation as generator of the time evolution of ip, see (3). On the 
other hand, the hermitian operator H, with real eigenvalues, was defined in a rather ad 
hoc manner, in order to absorb the non-hermitian terms found in (9). 

It turns out, however, that there is a straightforward, canonical, way to arrive at H. 
Indeed, the expectation value of H is simply the field Hamiltonian of the Dirac field, or 
in other words, the energy. 

In order to show this, we repeat some of the steps performed in 0. For convenience, 
we work in the Poincare gauge theory framework, i.e., the Lagrangian (1) depends on the 
fields ip, ip and their derivatives, as well as on the gravitational background fields T"■^^, e". 
Then, we consider the derivative of the Lagrangian density C = y/—gL, 

dip Olp^rn dip dlp^rn Dc^ 

Using the Dirac equations for ip and ip, we find 

[V^ t\h = -^.rr"L,. + T-e^,,, (13) 

where 

dL - BT 

dkiP + duiP - 5iL (14) 
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is the canonical stress-energy tensor of the Dirac field, while a^^™ = e'^SCj 5V°-'^^^ and 
T™ = —e~^5C/5e'^ are the spin density and the gravitational stress-energy tensor of the 
Dirac field (i.e., the sources of the gravitational field equations, see j2] for details). 
As outlined in [2], the field momentum is given by 



Vk= t\ dS,= t\ d'x, (15) 



and the integrated version of (13) reads 



in = - I Kr^^'k - TZel..k] d^^. (16) 



dt 

This is the usual form of a canonical conservation law, stating that Vk (for some component 
k) is conserved whenever the background fields and T'^''^ do not depend on the specific 
coordinate x^. (Quite similarly, if there were, in addition, a background Maxwell field, 
there would appear at the r.h.s. of (16) a term J \/—g j"^Am,k d^x, with the electric 
current density j^.) 

Especially, for the field Hamiltonian Ti = Vt = J y/—g t^t d^^; have 

^'H=- I V^g K.'-T^' - T^eU d\ (17) 



dt 

meaning that Ti is conserved whenever the background fields are time independent. These 
observations (and the fact that t\ can be generalized to include the gravitational fields 
themselves in a way that p*"*"' is conserved, see P]) are sufficient to identify Vk as the 
physical momentum of the system, and especially Ti as the field energy. 

We have already derived in the explicit expression for Ti in the fiat limit, it is the 
well known expression 

n = [ ip^HiP d^x = iiJ,HiP). (18) 



Therefore, in fiat space, the field Hamiltonian (which plays, in quantum field theory, the 
role of the generator of time translations), is simply given by the expectation value of 
the Dirac Hamiltonian H, i.e., of the generator of the time translations of the quantum 
mechanical wave function. 

From t*^ = (i/2)('^7*5tV' — dfipl^'ip), partially integrating using (6) and then the Dirac 
equation, or, alternatively, first the Dirac equation and then equation (9), it is not hard 
to show that in the general case, we find 

n = 1 ^^(7V)^V^ = {ip, Hi;), (19) 
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with H from (10). 

This result not only provides us with an interpretation of H , namely as the operator 
whose expectation values correspond to the energy of the system, but in addition, it 
justifies, a posteriori, the introduction of this arbitrarily defined operator. 

Consequently, the hermitian H should be interpreted as physical energy operator, 
while H has simply the mathematical interpretation of the time evolution operator. This 
interpretation can be further founded on the following. If we define, for some operator A, 
the operator dA/dt in the conventional way, i.e., by requiring 

(V',^^) = ^(^,^V^), (20) 
then, as is easy to show, with the inner product (7), we arrive at the usual equation 

-A^^^H,A]^^. (21) 

This relation, together with (3), clarifies the meaning of the operator H. 

The fact that both operators do not coincide is not really a problem, since in practice, 
the use of energy eigenstates makes only sense when the energy is conserved. This holds 
for time independent background fields, and in that case, we have H = H anyway. 

However, it can easily be shown that the momentum operator = id^ is not hermi- 
tian, but that, similar to (9), we have the relation 

= z[(p^^,^) - {iP^pM + (^, [a^ln ^J-g (7" ]^), (22) 
and consequently, in analogy to (10), the hermitian operator is given by 

i 
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Pt^ = Pt^+ In \/-gg^^- (23) 



Again, it turns out that the expectation value of this hermitian operator is identical to 
the field momentum from (15), i.e., 

V,= [V^ ^Hl'l')P,^ d'x = {iP,p,^), (24) 



in analogy to (19). In contrast to the situation for the Hamiltonian, the correction term 
in (23) will be different from zero even for the most simple field configurations. 

Finally, though unrelated to our discussion, we would like to remark that T"^ in (16) 
contains the full connection F^^j, and thus the full contortion tensor K"'^ (see for the 
explicit expression). Therefore, although only the axial part of K°-\ couples to the Dirac 
field in (1), the vector and tensor parts may nevertheless influence the field momentum 
through (16). 
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3 Non-unitary transformations 



In this section, we will make contact with the results previously obtained in literature. 
First, consider a transformation A acting on ip and on operators O as 

ij^ij' = Aij, O ^ O' ^ AOA-\ (25) 

where, at the same time, the scalar product (7) goes over to the flat scalar product, i.e., 

(V'i,V'2) = J ^lhV)^2 d'x ^ < ^^J,,^^J2 >= j d^x (26) 

We will refer to the tranformed quantities as being in the new representation. We require 
from A that it does not change physical quantities, i.e., eigenvalues and expectation values. 
Therefore, we must have 

(V'l,V'2)=<V'i,V'2>, (27) 

which leads to the condition 

A'fA = 7%*- (28) 
If the tetrad is in diagonal form, we have 7* = 607*^, and (28) this can be satisfied by 

A = {-9 9y^'. (29) 

If the tetrad is not diagonal, we first Lorentz rotate the vector e* such that = 0. 
(Then, of course, the gauge freedom is restricted, which is not unproblematic, because 
the Hamiltonian, as opposed to the Lagrangian, is not gauge invariant.) 

For A satisfying (28), all expectation values remain unchanged, especially the expec- 
tation value of H. Also, it is clear that if an operator O is hermitian with respect to ( , ), 
then, the transformed operator O' will be hermitian with respect to < , >. Therefore, 
(25) is a purely mathematical operation which does not affect the physics, just as are 
unitary transformations in special relativity. (In our case, the non-unitary part of A is 
absorbed by the simultaneous change of the scalar product.) 

The Dirac equation takes the form 

H'l/j' = AHA-^Ai) = Aidti) 

= idt^' -%{dtA)A-^ij\ (30) 

or, with the help of (29), 



idt^' = [H' + -dt In 9''W- (31) 
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If we compare this with (10), we see that in the new representation, the operator idt can 
be identified with the operator H, i.e., we have idt = H'. Clearly, just as H was hermitian 
w.r.t. the original scalar product, idt is hermitian w.r.t. the fiat scalar product. Thus, 
in the new representation, the energy operator is given directly by the time evolution 
operator. 

Note that the operator H', as before, is not hermitian, it satisfies equ. (9), where ( , ) 
is replaced by < , >. 

On the other hand, it has been shown that by the transformation 

H" = AHA-^ -iAdtA-\ = AtP, (32) 

where A still satisfies (28), the Hamiltonian H is brought into a form that is hermitian 
w.r.t. the fiat scalar product (see |Sllll)- It is not hard to show that this is indeed the case, 
but we wish to point out that there is more involved here then a mere transformation. 

First of all, it is clear that H" does not have the same eigenvalues as H. Thus, both 
operators are physically different. Secondly, as opposed to (25), which is applied to every 
operator, the transformation (32) acts specifically on the Hamiltonian. In the same way, 
you can bring, e.g., the momentum operator = into the hermitian form 

p';, = Ap^A-'-iAd^A-\ (33) 

which is clearly different from (32). Therefore, such transformations are not merely math- 
ematical tools, but are physical redefinitions of the operators H and p^. This becomes 
even more clear if we split the transformations into two steps. First, we perform the 
mathematical transformation (25), acting on all operators, but not changing any phys- 
ical quantities. Then, in a second step, we define H" = H' — iAdtA~^ (and similarly, 
p"^ = p' — iAd^A~^), which is simply the definition of new, hermitian operators. 

It is now clear that in the previous sections, we did exactly the same thing, namely 
we redefined H and in order to find hermitian operators (see (10) and (23)), with the 
only difference that we did it directly in the initial form, i.e., without performing the first 
step, which is completely unrelated to the hermiticity properties. 

To conclude, we showed that, apart from a purely mathematical change of represen- 
tation (25), whose purpose is to bring the scalar product into a fiat form (such that the 
hermiticity properties of the operators can be directly read off), the transformations (32) 
and (33) are exactly equivalent to the redefinitions (10) and (23). Therefore, the inter- 
pretation of (32) and (33) is that they are physically non-trivial redefinitions and not 
simply mathematical transformations. They are just as arbitrary as our equations (10) 
and (23), i.e., they can only be justified by physical arguments, especially by the fact that 
the resulting operators are hermitian. (Unfortunately, the transformation of the form (32) 
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is sometimes confused with a change of representation, and said not to change physical 
quantities, because it leaves the Dirac equation Hip = idfip form- invariant. However, this 
does not change the fact that H" is a different operator from H, as can be seen from 
the eigenvalues. This is because, with a change of the representation, the operator idt is 
transformed too.) 

The result of section 2, which relates the new operators to the field energy and the field 
momentum is therefore an a posteriori justification of the transformations (32) and (33) 
showing that the new, hermitian operators, really correspond to the physical quantities 
they are supposed to. 
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